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ABSTRACT 

We study the global efficiency of star formation in high resolution hydro dynamical 
simulations of gas discs embedded in isolated early- type and spiral galaxies. Despite 
using a universal local law to form stars in the simulations, we find that the early- 
type galaxies are offset from the spirals on the large-scale Kennicutt relation, and 
form stars two to five times less efficiently. This offset is in agreement with previous 
results on morphological quenching: gas discs are more stable against star formation 
when embedded in early-type galaxies due to the lower disc self-gravity and increased 
shear. As a result, these gas discs do not fragment into dense clumps and do not reach 
as high densities as in the spiral galaxies. Even if some molecular gas is present, the 
fraction of very dense gas (typically above 10^ cm~^) is significantly reduced, which 
explains the overall lower star formation efficiency. We also analyse a sample of local 
early- type and spiral galaxies, measuring their CO and Hi surface densities and their 
star formation rates as determined by their non-stellar 8/im emission. As predicted 
by the simulations, we find that the early- type galaxies are offset from the Kennicutt 
relation compared to the spirals, with a twice lower efficiency. Finally, we validate our 
approach by performing a direct comparison between models and observations. We run 
a simulation designed to mimic the stellar and gaseous properties of NGC524, a local 
lenticular galaxy, and find a gas disc structure and global star formation rate in good 
agreement with the observations. Morphological quenching thus seems to be a robust 
mechanism, and is also consistent with other observations of a reduced star formation 
efficiency in early- type galaxies in the COLD GASS survey. This lower efficiency of 
star formation is not enough to explain the formation of the whole Red Sequence, but 
can contribute to the reddening of some galaxies. 

Key words: galaxies: elliptical and lenticular, cD - galaxies: ISM - galaxies: star 
formation 



1 INTRODUCTION 

At least up to z ^ 1, galaxy colors follow a bimodal dis- 
tribution associated to their stellar morphology, with blue 
star-forming late-type spiral galaxies (LT Gs) and a red se- 
quence of early-type galaxi es (ETGs) fe.g. jKaulTmann et al.l 
l2003l : ICirasuolo et al.li20"07l ). Red ETGs harbor only very low 
levels of star formation, and alt hough these levels increase 
with redshift (|Tonini et al.ll201Ql ). ETGs remain well below 
the "Main Sequenc e " of s t ar formation (jNoeske et ahllioOTl; 



Salmi et 



Schiminovich et all l2007l : IWuvts et al] l201ll : 
2OI2I I. The conversion of cold gas reservoirs into stars on 



the scale of entire galaxies has long been considered to be 
a univers al process indepe ndent of the host galaxy proper- 
ties fe.g., iKennicutt lll998l V In that context, the quenching 
of star formation (SF) in ETGs is in general attributed to 
the quasi- complete removal and/or heating of cold gas reser- 
voirs, for instance because of merger-triggered starbursts. 
However, the residual gas content and late gas infall continue 
to fuel too high amounts of residual SF to popu late a red se- 
quen ce, unless extra mechanisms are invoked (|Gabor et al.l 

[mil)- 

Recent observations have actually shown that red ETGs 
can contain cold gas reservoirs. Massive and extended Hi 
reservoirs are often detect ed in or around ETGs th at are 
not in clusters (Morgant i et al ] I2OO6; Gros si et al ] [2009; 
pThom et al.ll2012|). In the volume-limited ATLAS^^ sample 
of nearby ETGs (|Cappellari et al.ll2011al , hereafter Paper I) , 
ISerra et all (|20ll hereafter Paper XIII) detect Hi in 40% 
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of the galaxies that are not in the Virgo cluster, with Hi 
masses from 10^ to a few 10^ Mq. More relevant for star 
formation, molecular gas is detected in 22% of ETGs, with 
an average molecular ^as-to-st ar mass fraction of a few per- 
cent fl-5%, lYoung et al]|201l l hereafter Paper IV), similar 
to what is found in some spirals. These molecular gas reser- 
voirs in ETGs are often found in relat ively compact discs 
([Crocker et al. ll201ll : [Davis et al.|[2011bl . hereafter Paper X), 
and as a result the typical surface densities of CO-traced 
molecular gas are not lower in ETGs than in LTGs. Even 
dense molecular trac ers, such as HCN , are now detected in 
red-sequence ETGs (jKrips et al.ll2010l : ICrocker et al.ll20ll 
hereafter Paper XI). The significant mass fraction and rela- 
tively high surface density of cold gas in many ETGs make 
their red optical colours more intriguing, as they should have 
specific rates of SF comparable to LTGs containing a few 
percent of gas. 

However, it has been recently recognized that SF may 
be a non-universal process that could depend on the prop- 
erties of the host galaxies. First, observations have sug- 
gested that the efficiency of star formation at fixed gas sur- 
face density incr eases from discs to mergers (Dad di et al.l 
I2OIQI : iGenzel et al.l |2Q10). Models explain this increase by 
the influence of galaxy-scale gravitati onal interactions on 
the sub-kpc-scale prop erties of the ISM (jTevssier et 111 120101 : 
'Bourn aud et al.ll201lh . Beyond the striking case of major 
mergers, there are also more systematic variations of the effi- 
ciency of SF and gas consumption timescale with host galaxy 
pro perties such as colors, stellar m ass or stellar mass density 
(see ISaintonge et~glll2011bl . I201I in the COLD-GASS sam- 
ple). The first studies of individual ETGs did not find evi- 
dence for departure from the usual Kennicutt relation estab- 
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lished for LTGs (jShapiro et al.ll201Ql : ICmcker et al. Il201lh . 

but there actually seems to be a clear tendency for lower SF 
efficiency and longer gas consumption timescale in galax- 
ies with redder colors, higher stella r mass concentrations , 
and / or higher stellar mass densities (ISaintonge et al.|[2011li , 
l2012l ). 

Theory has a promising explanation for the low SF ef- 
ficiency in ETGs that did not have their dense cold gas 
reservoirs entirely removed. The lower total disc self- gravity 
in the absence of a stellar disc, and the higher shear im- 
posed by stellar spheroids, can make a gaseous disc more 
stable in an ETG, provided that the gas velocity dispersions 
remain high enough. This increased disc stability was al- 
ready noted bv lKawata et al.l |2Q07l l for circumnuclear discs 
in elliptical galaxies (se e also early studies, for instance 
lOstriker Peebleslll973l . showing increased disc stability in 
presence of a spherical halo)^ 

It was shown by Martig et aP (|2009l , hereafter M09) in 
a cosmological zoom-in simulation that the morphological 
evolution of a galaxy from LTG to ETG can significantly 
reduce the star formation rate in the remaining gas reser- 
voir. This leads to red optical colors in a system that would 
remain blue if the star formation efficiency and gas con- 
sumption timescale were universal and independent on the 
host galaxy type. This mechanism, dubbed "morphological 
quenching" (MQ), has also been noted in zoom-in simula- 
tions at redshift 2; ~ 2, where the formation of big star- 
form i ng clumps is quenche d by bulge growth (jAgertz et al.l 
I2OO9I : ICeverino et al-llioTol ). 

In this paper, we present a series of idealised (isolated) 
high resolution (5 pc) hydrodynamic simulations of LTGs 
and ETGs (Section 2). In Section 3, we analyse the varia- 
tions of the SF efficiency of a given gas disc, depending on 
the host galaxy type. We confirm that on a Kennicutt dia- 
gram ETGs have a SF efficiency up to 5 times lower than 
LTGs at fixed gas surface density, consistent with a new 
analysis of the cosmological simulations from M09 (in M09 
the star formation efficiency of the simulated galaxy was not 
studied). In Section 4, we analyse observations of 8 nearby 
ETGs, for which we measure CO and H i surface density and 
star formation using their non-stellar 8/im emission. We ob- 
serve an effect quantitatively consistent with our theoretical 
predictions, namely a SF efficiency lower by a factor of ^ 2 
in ETGs with respect to a control sample of 12 spirals. The 
agreement between simulations and observations is further 
tested by a simulation aimed at mimicking the properties of 
NGC 524, a red ETG in our sample (Section 5). We discuss 
the robustness of the process and its potential cosmological 
implications for the growth of the Red Sequence in Section 6. 



2 SIMULATIONS 

2.1 Simulation in a cosmological context 

Morphological quenching is an effect that was first witnessed 
in a cosmological re-simulation presented in M09. In this 
Section, we briefly describe that simulation, but refer the 
reader to M09 for a detailed description of the simulation 
technique as well as the time evolution of stellar and gaseous 
properties for the simulated galaxy. 

We use a re-simulation technique that consists of first 
running a dark matter only cosmological simulation, and 



computing the merger and diffuse matter accretion history 
for one of the haloes. We then re-simulate this history at 
higher resolution using a Particle-Mesh code, in which gas 
dynamics is modelled with a sticky particle scheme. The 
spatial resolution is 130 pc, and the mass resolution is 
1.4x10^ Mq for stellar particles, 2.1x10^ Mq for gas parti- 
cles and 4.4x10^ Mq for dark matter particles. Star forma- 
tion is computed with a Schmidt law: the star formation rate 
is proportional to the gas density to the exponent 1.5, above 
a fixed threshold of 0.03 Mqpc"^. In this simulation, feed- 
back from supernovae explosions is not taken into account. 

The re-simulation is made by replacing each halo of the 
cosmological simulation by a new model galaxy made of a 
gas disc, a stellar disc and bulge and a dark matter halo, and 
by replacing each "diffuse" dark matter particle with a new 
group of gas and dark matter particles (mostly correspond- 
ing to accretion along filaments). The simulation starts at 
2; = 2 and follows the evolution of a galaxy down to 2; = 0, 
with mergers and gas accretion as prescribed by the initial 
cosmological simulation. 

The halo grows from a mass of 2x10^^ Mq at ^ = 2 to 
1.4x10^^ Mq at ^ = 0. During that time, the galaxy first 
undergoes a series of intense minor mergers that transform 
the initial z — 2 disc into a spheroid. From 2; :^ 1.2 down to 
z 0.2 is a much quieter phase, dominated by smooth dark 
matter and gas accretion, with an average gas accretion rate 
of 9 Mq yr~^. Finally, a major merger takes place at z = 0.2 
with a mass ratio of 1.5:1. 

During the quiet phase, a gas disc is gradually built 
within the ETG, with a cold gas mass growing from 
5.1x10^ Mq to 1.3 X 10^° Mq. In spite of this large gas 
content, the star formation rate remains at a low level, and 
the galaxy has a red color. This is the phase corresponding 
to a morphological quenching effect: the gas disc is stable 
against star formation because it is embedded in an ETG. 

In M09, we discussed the detailed properties of the gas 
disc, and analysed its stability. In this paper, we will show 
where the simulated galaxy stands on the Kennicutt relation, 
and will discuss the star formation efficiency as a function 
of host galaxy morphology. 



2.2 Idealized simulations of isolated galaxies 

2.2.1 General settings 

We perform a series of hydrodynamical simulations with the 
Adaptive Mesh Refinement (AMR) code RAMSES. The box 
size is 350 kpc and is covered with a 512^ grid, corresponding 
to a minimal level of refinement Imin = 9. Additional levels of 
refinement are added up to Imax = 16, so that the maximal 
spatial resolution is 5.3 pc. A cell is refined if it contains 
a gas mass greater than 3.7x10^ Mq or it contains more 
than 20 particles. This refinement strategy, together with the 
choice of the pseudo-cooling equation of state, ensures that 
the Jeans length is resolved by at least 4 cells at all levels 
of refinement so that artificial fragmentation is prevented 
(Truelove et al. 1997). Dark matter particles have a mass of 
7500 Mq and stellar particles have a mass of 2500 Mq. 

The gas is described with a barotropic cooling model (or 
"pseudo-cooling"), which has been s hown to produce a rea l- 
istic structure for the simulated ISM fB ournaud et al.|[20ld V 
Star formation is modelled with a Schmidt law: when the lo- 
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cal gas density exceeds a threshold of 30 cm~^, the local star 
formation rate is equal to e*pgas/tff, where tff is the free-fall 
time at the density pgas and e* is the efficiency. The efficiency 
is chosen to be 0.1% and is calibrated so that the sini ulated 
spirals follow the observed Kennicutt law (|Kennicutt | [T998V 
The combined choice of a relatively low threshold and a low 
efficiency is made to remove any possible bias due to some 
gas clouds being at the resolution limit, especially in the 
outer disc (where the overall density is lower so that cells 
are less refined). The size of these cells limits the hierarchy 
of structures that can collapse, and imposing a low threshold 
for star formation mimics the fact that some of these cells 
might have contained sub-clumps that we do not resolve, and 
should hav e formed stars at a low level. 

As in iBournaud et al supernova feedback is 

modelled in a kinetic form, where 20% of the energy of su- 
per novae is injected in the g as as a radial kick (the rest of 
the energy is assumed to be radiated away) within a bubble 
of radius 20 pc. We use the blast wave model described in 
iDubois k Tevssied (|2008h . 

2.2.2 Relaxation of initial conditions 

Initially, a galaxy model is not perfectly in equilibrium, but 
the gas disc needs to reach equilibrium before we start study- 
ing its properties. To allow a faster relaxation of initial con- 
ditions we increase the resolution progressively. We start the 
simulation without star formation, with an isothermal equa- 
tion of state at T = 10^ K and a maximal resolution of 43 pc 
corresponding to Imax = 13. When a steady state is reached, 
we turn on the pseudo-cooling, keeping Imax = 13. The next 
steps are to gradually increase the resolution up to Imax = 16 
so that the full resolution is reached, and to finally turn on 
star formation and feedback. After this, the simulations are 
each run for 100 Myr. 

This method ensures that a steady state is quickly 
reached, and the star formation rate is nearly constant for 
the main part of the simulation (Figure [1)). 

2.2.3 Galaxy parameters 

We simulate the isolated evolution of four different galaxies, 
with various Hubble types and gas fractions. Two of these 
galaxies are elliptical galaxies, with gas fractions respectively 
of 4.5 and 1.3% (the gas fraction is defined as the ratio of 
the total gas mass in the disk to the total stellar mass of 
the galaxy, and includes all gas phases). The two others are 
spirals, with also the same gas fractions of 4.5 and 1.3%. 

These gas fractions are chosen because they are repre- 
sentative of the gas content of local ETGs in the ATLAS^^ 
sample (Paper IV). In turn, the simulated spirals, which have 
a stellar mass and size similar to the Milky Way, have a much 
too low gas fraction to be representative of local spiral galax- 
ies, but serve as a comparison in our study of disc stability. 

All galaxies have the same dark matter halo, which has 
a total mass of 2.5x10^^ Mq (corresponding to 3.33 xlO'^ 
particles) , and is niodeled with a Burkert p rofile with a core 
radius of 14 kpc (jSalucci Burkertll2000l l. It is truncated 
at 130 kpc. In all cases, the stellar component has a total 
mass of 5.6x10^° Mq (corresponding to 2.24 xlO^ parti- 
cles). The elliptical galaxies are modeled with a Hernquist 
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Figure 1. Time evolution of the star formation rate for the four 
idealised high resolution simulations. This shows that the SFR is 
nearly constant as a function of time, a consequence of a careful 
relaxation of the initial conditions. The slow decrease with time 
of the SFR for the spiral with a 4.5% gas fraction is simply due to 
slow gas consumption. By contrast the elliptical with a 4.5% gas 
fraction first goes through a 75 Myr phase of slightly increased gas 
disk fragmentation and increased SFR, before reaching a steady 
state. 



profile with a characteristic radius of 1.5 kpc truncated at 
40 kpc. The spiral galaxies are made of a 5 x 10^ Mq bulge 
and a 5.1 xlO^° Mq disc so that B/D = 0.1. The bulge fol- 
lows a Plummer profile with a characteristic radius of 1 kpc 
truncated at 2 kpc and the disc has a Toomre profile with 
a characteristic radius of 3 kpc truncated at 20 kpc. The 
gas discs have a mass of 2.5 xlO^ and 7.5 xlO^ Mq (corre- 
sponding to gas fractions of 4.5 and 1.3% respectively). They 
follow an exponential profile with a characteristic radius of 
1 kpc truncated at 2 kpc. Their velocities are initialized in 
rotational equilibrium with the whole galaxy, but with no 
velocity dispersion. 

The gas disc parameters (mass and sizes) have been cho- 
sen to be within the range of observed values for local ETGs. 
The average extent of molecular gas in the ATLAS^^ sample 
i s 1 kpc, with only 15% of gas discs reaching a size of 2 kpc 
([Pavis et al.l 12012 . Paper XIV). However, this observed size 
only concerns molecular gas, which is the most centrally con- 
centrated component, while our simulated discs include all 
gas phases. In addition, the range of surface densities that 
the simulations explore (from ^ 10 to a few 100 Mq/pc^) is 
quite typical of the observed densities. 

A gas mass of 2.5x10^ Mq is however reached by very 
few local ETGs, and in that case our choice of disc pro- 
file corresponds to very high central surface densities (a few 
100 Mq/pc^). This represents an upper limit to the gas 
content of local ETGs, and the other simulated elliptical 
is probably more representative of the local population. By 
contrast, the spiral galaxies that we simulate are not rep- 
resentative at all of the local spirals (both in term of gas 
fraction and of gas disc radius, which are both too small), 
but will be useful to perform comparisons with the ellipti- 
cals. 
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Zgas [Mo pC-2] Zgas [M© pc-^] 

Figure 2. Simulated ETGs and spirals on the Kennicutt relation. Left : for a cosmological resimulation. Right : for idealised simulations 
with a maximal resolution of 5 pc. In both cases the shaded area is the Icr zone from Kennicutt (1998) and the solid lines are fits to 
the simulated data points (for the idealised ellipticals, a fit is performed separately for the two galaxies since they have very different 
behaviours). The dashed red lines are there to guide the eye and are simply the fits for the spiral population but with an efficiency divided 
by 2.5 (left) or 5 (right). 



3 GAS DYNAMICS AND STAR FORMATION 
IN SIMULATED GALAXIES 

3.1 Star formation efficiency 

Hereafter, we define the star formation efficiency a as : 

5]SFR = a X Sfas , (1) 

where the star formation surface density Ssfr is expressed 
in units of Mq yr~^ kpc~^ and the gas surface density Sgas 
is in Mq pc~^. This definition of the efficiency is chosen 
because it is the easiest one when the exponent (3 changes 
from one galaxy to another. 

On the left panel of Figure [2l we show the gas and star 
formation rate surface densities in the cosmological simu- 
lation studied in M09. We measure these quantities during 
the morphological quenching phase in 5 snapshots (500 Myr 
apart) of the simulation. For each snapshot, average values 
of Sgas and XlsFR are computed within the optical radius of 
the main galaxy (i?25 in g band), and in 10 annuli centered 
on the center of the galaxy (each annulus has a width of 
1 kpc, and the galaxy is taken face-on). The same method 
is applied to other snapshots of the same simulation, when 
the galaxy is actually a star forming disc at 2; ^ 2 before the 
first minor mergers occur, and at ^ ^ 0. 

Except at very low gas density, the star formin g discs 
lie within la of the best fit proposed bv iKennicutt I jl998). 
The red ETGs undergoing morphological quenching roughly 
occupy the same region of the Kennicutt plot as the star 
forming discs, although very high gas densities are not found 
in the MQ phase. As expected, on average the quenched gas 
discs tend to lie slightly below the star forming galaxies (i.e., 
to have a lower User for a given Sgas). We fit the efficiency 
and slope of the Kennicutt relation (Eq. [!)) for each series 
of points independently, not taking into account the points 
at low gas density and low star formation rate that would 
artificially skew the slope to a too high value (and might be 
affected by our choice of threshold for forming stars in the 
simulation) . We list in Table [1] the best fit values for the 
efficiency and slope of the Kennicutt relation. The efficiency 



Table 1. Best fits values for the efficiency and slope of the Ken- 
nicutt relation for simulated galaxies 



Simulation 




Efficiency 


Slope 


Cosmological 


Disc 




1.9 X IQ-^ 


1.4 


Early-type 




6.5 X 10-6 


2.3 


Idealized 


Spirals 




1.7 X 10-4 


1.6 


Elliptical, f^= 


=4.5% 


1.4 X 10-4 


1.6 


Elliptical, f^= 


=1.3% 


7.9 X 10-6 


2.0 



and slope for the spiral galaxies are ver y close to the best 
fit values proposed by IKennicutt I (|l998l ). while the ETG 
has much lower efficiency but also a much steeper slope. 
This steep slope makes it hard to compare the efficiency 
directly with the one found for the spirals. To help with the 
interpretation of the results, we draw a line on Figure [2] with 
the same slope as the spirals' one and that goes through 
the points corresponding to global averages. This line gives 
an efficiency of star formation 2.5 times lower in the ETG 
compared to the spirals. 

We test if this result still holds in the 5-pc resolution 
simulations of gas discs embedded in an idealised spiral or 
elliptical galaxy. We now measure the gas and SFR surface 
densities within 500 pc wide annuli between and 2 kpc from 
the center of the galaxy (see right panel of Figure [2]) . This 
amounts to four data points per simulated galaxy, except 
for the spiral with the lowest gas mass: in that galaxy the 
outer regions have a very low gas density and do not host 
any star formation, so that only the first two radial bins are 
considered. 

We fit both simulated spirals with the same slope and ef- 
ficiency. We perform separate fits for the two elliptical galax- 
ies since they have a very different behaviour (see the result- 
ing best fit values in Table [1]). The elliptical with a large gas 
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Table 2. Gas disc properties and star formation rate (averaged over the last 5 Myr of the simulations) for the four idealised simulations. 



Simulation 


(7± [km/s] 


[km/s] 


ig{> 10^ cm-3) 


f^(> 10^ cm-3) 


Mg (> 10^ cm-3) [M©] 


SFR [Mo yr-i] 


Elliptical - f^=1.3% 


0.02 


0.04 


0.26 


0.00 


0.0 


0.1 


Spiral - f^=1.3% 


2.48 


5.22 


0.89 


0.51 


3.8 X 10^ 


3.8 


Elliptical - f^=4.5% 


1.05 


1.74 


0.72 


0.18 


4.5 X 10^ 


4.7 


Spiral - f^=4.5% 


2.97 


5.26 


0.87 


0.43 


10.7 X 10^ 


11.3 



fraction (and thus a high gas surface density) lies very close 
to the spirals' best fit : the slope is the same in both cases, 
and the star formation efficiency is only 1.2 times lower. By 
contrast, the elliptical galaxy with a lower gas content is 
clearly offset from the relation followed by the spirals, and 
has a lower star formation efficiency. As for the elliptical in 
the cosmological simulation, the slope of the Kennicutt re- 
lation is also steeper for this galaxy. Drawing a line passing 
through the data points but with the same slope as for the 
spirals gives a factor of 5 difference between the star forma- 
tion efficiency in this elliptical compared to the spirals. 

We thus do not predict that ETGs are completely off 
the relation observed for spirals, but that they have a star 
formation efficiency up to a factor of 5 lower than spirals. The 
amount by which the efficiency is reduced depends strongly 
on the gas disc properties, which we will study next. 

3.2 Gas disc properties 

To better understand the difference in star formation effi- 
ciency as a function of the host galaxy morphology, we study 
the properties of the gas discs in the four generic simulations. 
The top panels in Figure [3] show the gas surface density maps 
after 100 Myr of evolution at full resolution. We also show 
the gas density profiles and the gas density probability dis- 
tribution functions (PDFs) in Figure [H While all simulations 
started with a smooth exponential gas disc, we find that the 
subsequent evolution depends strongly both on the total gas 
content (more precisely on the gas surface density) and on 
the Hubble type of the host galaxy, in agreement with our 
previous results on morphological quenching. 

For a relatively low surface density and gas fraction, we 
find that the gas disc is extremely stable when embedded 
in the elliptical galaxy. It keeps an exponential profile with 
a similar scale-length as the initial disc, and does not frag- 
ment into dense clumps, although some small spiral armlets 
can be seen in Figure (3] The gas densities do not reach val- 
ues greater than 10^ cm~^, and only 26% of the gas mass is 
found at densities greater than 10^ cm~^ (a typical value for 
molecular gas formation). As a consequence, the star forma- 
tion rate remains very low, at only 0.1 Mq yr~^. Note that 
the gas disc is stable in spite of a very low turbulent velocity 
dispersion (Tableland bottom left panel of Figure [3]). The 
low velocity dispersion is the consequence of both the gravi- 
tational stability and very few supernovae explosions. Some 
of the gas disc stability is probably provided by the thermal 
velocity dispersion. 

By contrast, we find that a spiral galaxy with a similar 
initial gas content has very different behaviour. The gas disc 
is heavily fragmented, and the migration of the gas clumps 
under the effect of dynamical friction concentrates the gas 



towards the center of the galaxy so that the final density 
profile is strongly changed (Figures [3] and 3]). Within these 
clumps, the gas can reach densities of a few 10^ cm~^ (see 
the PDF in Figured]). Contrary to the case of the elliptical 
galaxy, 89% of the gas mass is found at densities greater than 
10^ cm~^, and 51% of the mass reaches densities typical of 
prestellar cores (10^ cm~^). The spiral galaxy then forms 
stars at an average rate of 3.8 Mq yr~^, a rate nearly a 
factor of 40 greater than in the elliptical galaxy. 

The differences between spiral and elliptical are not so 
pronounced at high gas surface densities. We find that in 
the elliptical with a 4.5% gas fraction, the disc is not stable 
against fragmentation into clumps (see Figure [3]). Both the 
surface density map and profile are rather similar for the 
spiral and the elliptical. However, some differences can still 
be seen in the gas density PDFs, with a PDF shifted to higher 
densities in the spiral galaxy. Indeed, while the spiral and 
the elliptical have a similar fraction of their gas in molecular 
form (72% in the elliptical 87% in the spiral), the fraction of 
very dense (> 10"^ cm~^) gas is 2.4 times higher in the spiral, 
resulting in a SFR 2.4 times higher than in the elliptical. 

In the idealised simulations the SFR is actually directly 
proportional to the mass of gas at densities greater than 
10^ cm~^ (see the last two columns of Table[2]), except in the 
elliptical with the smallest gas fraction. In that case there 
is some residual star formation (at 0.1 Mq yr~^) in spite 
of no gas at high densities; this is because of the very low 
density threshold (30 cm~^) we used for star formation in 
all simulations. 

The effect of morphological quenching is thus to change 
the shape of the gas density PDF in ETGs, and to change 
their fraction of dense gas. Given that we assume a star 
formation law which is universal on local scales, the smaller 
amounts of dense gas naturally translate into lower levels 
of star formation in ETGs. We will discuss this further in 
Section 16.21 but we will first study if such a morphological 
quenching effect can also be observed in local ETGs. 



4 RESOLVED OBSERVATIONS OF 
MOLECULAR GAS AND STAR 
FORMATION IN EARLY-TYPE GALAXIES 

To compare our simulations to observed galaxies we select 
eight galaxies from the ATLAS^^ sample of early-type galax- 
ies. All eight have molecular gas maps based on interferomet- 
ric CO (1-0) data, atomic gas maps, and star formation rate 
maps obtained using the 8/im Polycyclic Aromatic Hydro- 
carbon (PAH) band as a star formation rate indicator. This 
sample is made of all galaxies for which this information 
is available, after removing one galaxy (NGC 3489). NGC 
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Figure 3. Face-on maps of the gas surface density and velocity dispersion. The top panels show the gas surface density (5x5 kpc 
panels, colormap in units Mq/pc^) for the 4 generic simulations. The bottom panels show the corresponding in-plane velocity dispersion 
cr|| = Y^cr| + c''^ 5 mass-averaged along the line of sight, computed in 80x80 pc pixels (5x5 kpc panels, colormap in units km/s) 



c77l03 



10^ 



! 10^ 

OS 



10° 



10° 







— Initial profile 






— Spiral - fgas=1.3% 






— Elliptical - fgas=1.3% ' 









1 2 3 

Radius [kpc] 



— . — — . — — . — — . — — . — , 


— . — — . — — . — — . — ™. 


Initial disk 




Spiral - fgas=1.3% 




Elliptical - fgas=1.3% 






^^^^^^ 


. , 1 . ... 1 . ... 1 . .'■ . , 1 . ... 



10° 



10° 



— Initial profile 
Spiral - fgas=4.5% 

— Elliptical - fgas=4.5% 




3 


1 2 3 
Radius [kpc] 




4 E 




T • ^^'T • ■ ■ ■ 1 ^' ^^'T • ^^'T • — 

/'■\ - 


Initial disk 






/ \ 


Spiral - fgas=4.5% 




1 1 

/ / 
1 / 


Elliptical - 


fgas=4.5% 




/ / 
/ / 

/ /\ / 1 / 

/ 1 / 
/ / 1 / 








/' 1 








1 / ,'..1- ..1 . . .,1 . 







10-2 10-1 10° 10^ 102 10^ 10"^ 10^ 10^ 10^ 
Gas density [cm"^] 



10-2 10-1 10° 10^ 10^ 10^ 10"^ 10^ 10^ 10^ 
Gas density [cm"^] 



Figure 4. Gas surface density profiles (top panels) and gas density probability distribution functions (PDF, bottom panels) for the spiral 
and elliptical with a gas fraction of 1.3% (left panels) and 4.5% (right panels). In each case the initial profile is also shown. 
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Table 3. The sample of observed ETGs 



Name 


Distance 


H2 mass 


H I mass 




PA 


Axis Ratio 




(Mpc) 


lofffMrT^) 


loff(Mr7^) 

oV / 










Lo.o 


7 S 
< .0 






OO.O 


±.uo 


NGC 2768 


21.8 


7.8 


7.8 


0.258 


180 


1.20 


NGC 3032 


21.4 


8.7 


7.8 


0.293 


92.5 


2.06 


NGC 4150 


13.4 


7.7 


6.0 


0.264 


146 


1.44 


NGC 4459 


16.1 


8.2 


< 6.6 


0.264 


89 


1.47 


NGC 4477 


16.5 


7.4 


< 7.0 


0.261 


227 


1.20 


NGC 4526 


16.4 


8.8 


< 7.9 


0.268 


108 


2.28 


NGC 4550 


15.5 


6.9 


< 6.5 


0.268 


179 


3.55 



Distances f rom Paper I. H 2 mass es from lCrocker et al. ] (|201ll ) and use Xco = 3 x 10^*-* cm ^ (K km s ^. masses from 
Paper XIII or I Crocker et al.1 (|201ll) . /a.e is the multiplica tive factor applied to the IRAC SS/j^m image when it is subtrac t ed ofl P the 
8/im image in order to remove the stellar emission at 8fim (|Shapiro et al.|[201o|). Position angles are from ICappellari et al ] (l2007h and 
axis ratios are derived from best inclinations from bavis et ahl (|2011al . Paper V). 



3489 has a very unusual molecular gas distribution, includ- 
ing spiral arms. This suggests that the gas is gravitation- 
ally unstable, and is consistent with NGC 3489 being a flat, 
discy galaxy, with stars that might contribute to the disc 
instability. In addition, the molecular gas has an irregular 
velocity field su^ ^esting that it h as been recently accreted 
or disturbed (Cro cker et al. I [201 1'). Because of these uncer- 
tainties, we exclude this galaxy from the sample. Details of 
the sample galaxies are listed in Table O 

For comparison, we analyse a sample of twelve late- 
type galaxies: NGC 628, NGC 2903, NGC 3351, NGC 3521, 
NGC 3627, NGC 4736, NGC 4826, NGC 5055, NGC 5194, 
NGC 5457, NGC 6946, and NGC 7331. This sample con- 
sists of all late-type galaxies with available CO from BIMA- 
SONG (selecting only tho se with both BIMA and 12m obser- 
vatio ns lHelfer et al.ll2003l ). Hi fro m THIN GS (Walter et al] 
I2OO8I ) and Sum data fr om SINGS ( Kennicutt et al.ll200 3) or 
LVL (jPale et all [200^1. These galaxies are at similar dis- 
tances (5-15 Mpc) to our ETG sample galaxies. 



4.1 Observational gas surface densities 

Gas surface densities are derived using CO and Hi inter- 
ferometric observations. We use the CO (1-0) transition 
as a tracer of molecular gas, using a conversion factor of 
Xco = 3 X 10^° cm-2 (K km s-^)-^ Half of the ETG 
CO maps (NGC 3032, NGC 4150, NGC 4459, NGC 4526) 
were obtained using the Berkeley Illinois Maryland Array 
(BIMA ) with details of ob servational setup and data reduc- 
tion in I Young et~aD (|2008l ). The other half of the CO maps 
(NGC 524, NGC 2768, NGC 4477 and NGC 4550) were 
obtained with the Plateau de Bure I nterferometer (PdBI) 
wit h details in ICrocker et al.1 (|2008l ). ICrocker et aD (2009) 
and lCrocker et al. I 1^201^ . 

Our choice of Xco is based upon the typical value 
for Milky Way clouds. It is known that Xco varies based 
upon metallicity, reaching much high er values in low- 
metallicity systems (e.g. Wilson 1995; Ar imoto et al. |[l996V 
Meanwhile, highly active galaxies such as starbursts and 
ULIRGS have lower values of Xco (e.g. IWild et al.1 Il992l : 
ISolomon etHI Il997l l. Unfortunately, little is yet known 
about what the appropriate value for early- type galaxies may 
be. Recent theoretical modeling suggests that Xco is not a 



Table 4. Observational results for early- type galaxies 



Radius 


^gas 


^SFR 


kpc 


M0 pc-2 


Mq yr~^ kpc~^ 


NGC524 






0.31 


77.9 ±16.3 


0.0373 ± 0.0051 


0.62 


32.8 ±6.6 


0.0195 ± 0.0024 


0.92 


11.4 ±2.3 


0.0084 ± 0.0013 


NGC2768 






0.18 


47.4 ±9.5 


0.0173 ± 0.0038 


NGC3032 






0.58 


204.1 ±39.6 


0.1029 ±0.0032 


1.15 


81.3 ±15.1 


0.0218 ±0.0007 


1.73 


20.4 ±3.5 


0.0037 ±0.0001 


NGC4150 






0.44 


77.1 ± 15.7 


0.0380 ±0.0017 


0.89 


7.2 ±1.5 


0.0035 ± 0.0002 


NGC4459 






0.57 


119.1 ±23.9 


0.0632 ± 0.0034 


1.14 


17.3 ±3.5 


0.0051 ±0.0007 


NGC4477 






0.22 


114.5 ±23.0 


0.0141 ±0.0045 


NGC4526 






0.35 


637.4 ± 127.6 


0.2037 ± 0.0095 


0.70 


265.3 ±53.1 


0.0705 ± 0.0033 


1.05 


45.7 ±9.2 


0.0094 ± 0.0012 


NGC4550 






0.33 


15.5 ±3.2 


0.0081 ±0.0017 



function of galaxy morphology, depending most strongly on 
metallic ity and more weakly upon radiation field and column 
density (jNaravanan et al.ll2Q12 : lFeldmann et al.ll2012 l. None 
of the early-type galaxies chosen for this work are starbursts 
and furthermore, many have low [Oiii]/H/3 ratios signalling 
relatively high gas-phase metallicities. Therefore a choice of 
^co = 3 X 10^° cm~^ (K km s~^)~^ is at least plausible for 
the early-type galaxy sample and simplifies the comparison 
with the spiral galaxies, for which we will use the same value 
of Xco • 

We sum the fiux density in the CO maps of these galax- 
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ies within annular ellipses. The annuli are centered at each 
galaxy's optical center with widths equal to the average of 
the interferometric beam major and minor axes in order to 
keep the measurements independent. The position angle and 
axis ratio of the ellipses are listed in Table [S] No background 
subtraction is performed for the CO photometry. 

The Hi data obtained with the Westerbork Synthe- 
sis Radio Telescope are less well resolved (beams typically 
> 30^^) and so we only have 'central' values of H i mass (Pa- 
per XIII). We assume this mass is uniformly distributed over 
the observed CO distribution. While the H1/H2 ratio almost 
certainly varies within each galaxy, we note that the inclu- 
sion of the H I is always a minor consideration for the early- 
type galaxies, as they have only ^ 10% or less of their cold 
gas mass in H i in these regions. 

An identical approach is followed for the late- type galax- 
ies, including using the same Xco factor. The only change 
is the use of the interferometric Hi maps from the VLA as 
obtained by the THINGS project. The resolution of these 
maps is ~ 6^^ so we are able to use the Hi map itself, com- 
bined with the CO maps to obtain a total cold gas surface 
density. 

The dominant error in the gas surface density is the 
calibration uncertainty associated with the CO observations, 
taken to be 20%. We also add in quadrature a component 
for the estimated photometric error, based upon the noise 
level in the original data cubes. 

4.2 Observational star formation rate surface 
densities 

For both early and late- type galaxies, we compute star for- 
mation rates using their non-stellar 8/im emission. Emission 
from this wavelength predominantly traces band emission 
from PAHs. We choose S/xm as a star formation rate indi- 
cator because its resolution is better than some other avail- 
able indicators (2.82'' compared to: GALEX FUV at 4.5'', 
Spitzer 24/im at 6.4"). We note that Ha would have a better 
resolution than 8/im, but too few (2) of our ETGs have avail- 
able Ha maps. To obtain a measure of the non-stellar 8//m 
emission, we subtract a fraction of the 3.6/im IRAC band 
(/s.e; dominated by stellar continuum) from the 8/im IRAC 
band. The fraction used changes from galaxy to galax y and 
i s bas ed on stellar population models as in Shapiro et al.l 
(|201QI 1. These values are listed for the early- types in Table [3] 
and a constant fraction of 0.293 is used for the late- type 
galaxies. 

Photometry is performed in the the same elliptical an- 
nuli as used for the gas. Background subtraction is performed 
before the stellar subtraction in both the 3.6 and 8/im images 
and the recommended aperture corrections are applied. 

For t he conversion to SFR, we turn to the papers of 
IWu et al.1 (|2005al ) and IZhu et al.1 (|2008l ) which give SFR 
conversions for the IRAC 8/xm band based upon compar- 
isons to other SFR tracers (radio. Ha, IR, UV). We con- 
vert these from a Salpeter to Kroupa Initial Mass Func- 
tion (IMF) and find conversion factors from 1.06 xlO~^^ 
to 1.303 xlO~^^ (converting from luminosity in erg s~^ to 
SFR in Mq yr~^), based on the tracers used for compar- 
ison. We opt to take a roughly average conversion factor 
of 1.2xl0~^^ in this work. We note that some fraction of 
the 8/im emission will be contributed by old stars exciting 
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Figure 5. Star formation rate and total gas surface densities for 
a sample of local late-type (blue squares) and early-type galaxies 
(red dots). For each galaxy, these values are computed in several 
radial bins depending on the extent of the molecular gas. We also 
show the Icr zones from Kennicutt (1998) (shaded area) and Bigiel 
et al. (2008) (dotted lines), both corrected for our Xqq value and 
for a Kroupa IMF. The solid blue and red lines are fits to the 
data points, showing an average star formation efficiency lower 
by a factor 2 in ETGs. 

PAH molecules and n ot exclusively UV ph otons from young 
star- forming regions (jCrocker et al.|[2013l ). Due to the con- 
centration of old stars in the bulges of ETGs (where the 
molecular gas is also found), the early- type galaxies likely 
have a larger c ontribution to their 8/im emission fr om old 
star heating ((Xilouris et al.ll2004l : [Kaneda et al.ll2008l ). How- 
ever, the 8/im emission is w ell tied to the hot du st emission 
in the 24/im band in ETGs (|Crocker et al. Il201lh . Thus the 
8/im SFRs are not likely to be far off, but should be strictly 
regarded as upper limits for the early- type galaxies. 

The error for the derived star formation rates includes 
the photometric error, a 3% error from IRAC calibration in 
both the 3.6 and 8/im images and a systematic error of 10% 
on the stellar subtraction fraction. These are combined in 
quadrature. 

4.3 Results 

Our derived gas and star- format ion surface densities are 
listed in Table S] for the early-type galaxies. In Figure O 
we plot the values for both the early- and late- type galaxies. 
While there is some spread in both populations, the points 
corresponding to the annular averages for early-type galaxies 
are biased towards having lower values of Ssfr at a given 
gas surface densi ty. We note that this is an improvement 
on earlier works (|Shapiro et al.l I2OIOI : I Crocker et al. I l201lh 
that were unable to distinguish an ETC offset. We now spa- 
tially resolve the gas and star formation and we also gain 
diagnostic power by directly treating spiral galaxies with an 
identical methodology. Additionally, as noted above, we ex- 
cluded a galaxy with an irregular molecular gas distribution, 
which the previous studies had included. 

We fit each population separately, and find that the fits 
for ETGs and LTGs show a similar slope, whic h is close to 
linear as previously found for LTG samples fe.g. jBigiel et al] 
I2OO8I ). However, the efficiency of star formation is lower by 
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a factor of ^ 2 in the ETGs, and all ETG data points are 
below the best fit line for the LTGs. This lower efficiency of 
star formation is consistent with morphological quenching 
and the findings from the simulations. 



5 A DETAILED STUDY OF NGC 524 

In our series of idealized simulations, we studied gas disks in 
perfectly spherical, non-rotating elliptical galaxies. However, 
in the local Universe such galaxies are very rare, and rarely 
contain molecular gas (Paper IV): most ETGs containing 
molecular gas are fast rotators. To understand if morpho- 
logical quenching is still efficient in a more realistic galaxy 
with a stellar disk component, and to test the agreement be- 
tween observations and simulations, we have also performed 
a RAMSES simulation dedicated to the study of one of our 
sample ETGs: NGC 524. 

NGC 524 is a red lenticular galaxy (a fast rotator) seen 
close to face-on, it has a central disc containi ng 6.7x10^ M^t^ 
of m olecular gas within a radius of 1.1 kpc (jCrocker et al. I 
[2011), and an estimated SFR of 0.035 - 0.051 Mq yr"^ 
NGC 524 was chosen for this comparison because of the reg- 
ular, symmetric structure of its gaseous and stellar compo- 
nents (with also regular kinematics), suggesting the absence 
of recent accretion events. It also hosts a well-resolved dust 
disk, the structure of which we can compare to simulations. 



5.1 Simulation setup 

The initial conditions for the simulation of NGC 524 
have been derived based on an individ ual Mult i- Gaussian- 
Expansion model (Emsellem et al.'"l994") in the course of the 
ATLAS^^ project (Scott et al. 2012). 

The MGE model is fitted to gro und-based i mages with 
the IDL fitting implementation from lCappellari I ([2002 ) : the 
stellar component is defined as a sum of two-dimensional 
Gaussians which are then deprojected to three-dimensional 
Gaussians using an inclination of 19° for NGC 524. We find a 
total stellar mass of 4.3 x IO^^Mq. The gas and dark matter 
components are finally added using an MGE representation 
of an exponential disc and a Burkert profile, respectively. 
The dark matter halo has a core radius of 15 kpc (corre- 
sponding to 3i?e), and its mass is 1.6 x IO^^Mq so that dark 
matter represents 1/3 of the total mass within 3i?e- 

The full model (stars, gas, dark matter) is then used 
to generate initial conditions using the python pymge mod- 
ule (Emsellem et al., in prep). The pymge package takes the 
MGE models and solves Jeans Equations directly, assum- 
ing a cylindrically oriented velocity ellipsoid, while allowing 
for general axis ratios of the ellipsoid and concatenation of 
gaussians co mponents. This i s done with the efficient MGE 
formalism of lEmsellem et al. (Il994[ ), with the general axis 
ratio case as in ICappellaril ( 2008l l, and grouping of Gaus- 
sians as described in Emsellem et al. (in prep). We thus 
derive a model including the positions and velocities for 2.1 
million particles (10^ for the stars, 10^ for the halo, and 10^ 
for the gas) consistent with the specified underlying MGE 
mass distribution. Note that the gas component is replaced 
when feeding these models in the RAMSES code to make 
it consistent with the AMR scheme. However, the gas itself 
is included during the generation of the initial conditions to 



make sure we obtain a fully self- consistent dynamical reali- 
sation. 

In the RAMSES simulation, the gas disc initially has a 
mass of 1.6 X 10^ Mq , with an exponential scale- length of 480 
pc. It is initially truncated at 2.4 kpc. For this simulation, 
the refinement strategy, cooling model, star formation and 
feedback parameters are the same as for the generic idealised 
simulations. 



5.2 A comparison of gas disc properties and star 
formation rate 

We show in Figure [6] a comparison of the simulated gas 
disc morphology with a map of the CO emission obtained 
with the Plateau de Bure Interferometer and with a F555W 
unsharp-masked image from the Hubble Space Telescope 
Wide-Field Planetary Camera 2, showing tightly wound spi- 
rals of optically obscuring dust. This tightly wound spiral 
structure is also found in the simulated gas disc, which has 
a fiocculent structure with many small spiral armlets. The 
small pitch angle is consistent with the gala xy having a mas- 
sive bulge, and thus a high shear rate ( Julian Toomra 
ll966l : lGrand et al.ll2012l ). 

The total star formation rate is also similar in both 
cases, with a simulated value of ^0.045 Mq yr~^, while the 
observed SFR is in the range 0.035 - 0.051 Mq yr-^ This 
low value is consistent with the simulated gas density PDF, 
which is typical of a stable gas disc, with a maximal density 
of 5x10^ cm-^ 

The star formation and gas surface densities are com- 
puted for the simulation in the same radial bins as for 
the observation, and the result is shown in Figure [3 We 
see an offset between observation and simulation on this 
Kennicutt plot, mostly because the simulation does not re- 
produce the observed gas density profile, which is much 
flatter. However, both highest density points are a good 
match. Note that some differences between model and ob- 
servation may also arise because NGC 524 is seen close to 
face-on, which makes t he mass deprojection very degenerate 
(iLabla nche et al.| |2012|. Paper XII). 

However, no simulation parameter was tuned to repro- 
duce the observed SFR: the parameters are kept the same 
as for the generic simulations, including the star formation 
efficiency (that was initially adjusted to reproduce the po- 
sition of disc galaxies on the Kennicutt relation). The good 
match we obtain is reassuring both from the numerical and 
observational sides. 



6 DISCUSSION 

6.1 Robustness of the morphological quenching 
mechanism 

We witnessed and studied morphological quenching in dif- 
ferent types of simulations, from low-resolution cosmologi- 
cal zooms to high-resolution idealised galaxies. In addition 
to having different resolutions, these simulations were per- 
formed with different gas modelling techniques and different 
implementations of sub-grid recipes for star formation and 
feedback. We also note that while the idealised case assumes 
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Figure 6. Gas distribution in NGC 524. Left : simulated face-on 
gas surface density map (4x4 kpc). Right : integrated inten- 
sity map of CO(l-O) (Jy beam"-*^ km s"-*^), revealing the distri- 
bution of cold gas in NGC 524 (top) and unsharp-masked HST 
WFPC2 F555W image showing optically obscuring dust, with 
CO (1-0) contours overlaid in black (bottom). A similar morphol- 
ogy is found in the model and the data, with a tightly wound 
spiral structure in both cases. 
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Figure 7. A comparison of observed and simulated star formation 
and gas surface densities for NGC 524, computed in three radial 
bins with outer radii of 0.3, 0.6 and 0.9 kpc. The solid line repre- 
sents Kennicutt (1998) best fit, and the dashed lines show the 1-cr 
range. The gas density profiles do not match, with a much steeper 
profile in the simulated galaxy. However, the overall decrease of 
the star formation efficiency in the simulation is consistent with 
the observed value, and the total SFRs are very close. 



a perfectly spherical galaxy, both the cosmological simula- 
tion and the simulation of NGC 524 (a lenticular galaxy) 
reproduced more realistic mass distributions. MQ is thus ef- 
fective as long as galaxies are bulge- dominated, and not only 
in perfectly spherical galaxies. 

Finally, MQ has also been observed in simula tions 
by other groups, in p articular by 'Agertz et ajj (|2009l l and 
ICeverino et al.l (|201Q| 1 f or hi^h r edshift galaxies (see also a 
theoretical discussion in lPekel et al.l t2009V This means that 
from a simulation point of view, this is a robust mechanism, 
with now several different simulations showing a reduced star 
formation efficiency in ETGs. 

From the observational point of view, several factors 



could influence the position of our sample on the Kenni- 
cutt relation. The offset between ETGs and spirals could in 
theory arise from a different conversion of CO luminosity 
into mass of molecular gas for these two galaxy types. How- 
ever, removing the offset would require the ETGs to have 
a lower Xco than spirals, which is inconsistent with their 
higher metallicities. There is also some uncertainty on the 
derivation of the star formation rate from the non-stellar 
8/im emission, but the choices we made should be conserva- 
tive (see Section 4), and provide upper limits to the positions 
of our ETGs on the Kennicutt relation. 

Some final uncertainties concern the IMF in the 
quenched gas discs in the ETGs. Previous studies found 
that the ongoing star formation events in more mas- 
sive and denser galaxies have a more top-heavy IMF 
than Salpeter. They used indirect methods relating galaxy 
colour or UV fluxes to the Ha ionized ^as emis- 



sion ([Hoversten Glazebroo^ l2008l : iMeurer et all l2009l; 
iGunawardhana et alj 201 ih . More direct indication of the 
non-universality of the IMF in galaxies comes from the 
analysis of IMF-sensitive spectral absorption features. 
These indicate a bottom -heavy IMF in massive ellipticals 
([van Dokkum Conrovl[2010 V The dynamical modelling of 
the ATLAS^^ sample found a systematic variation in the 
mass normalization of the IMF within the ETG population. 
This varies from Kroupa/Chabrier to Salpeter or heavier, 
as a function of the ga laxy density or velocity dispersion 
(ICappellari et alJ [ioTjal lbl, hereafter Paper XX). The spec- 
tral and dynamical analyses measure the integrated IMF 
over the past history of the galaxy and not the present day 
IMF that is relevant for this study. However all these results 
together indicate that the IMF varies systematically with 
galaxy properties and possibly with time. If less high-mass 
stars, or more low-mass ones, are formed in the present-day 
ETGs with respect to spirals, the same luminosity would 
correspond to a higher mass of stars formed, and we would 
thus underestimate the SFR in ETGs (for instance by a fac- 
tor of 1.44 if ETGs had a present-day IMF closer to Salpeter 
than to Kroupa). 

The relatively good agreement between simulation and 
observation of NGC 524 gives us confidence that the reduced 
SF efficiency in ETGs is real and not an artifact of the han- 
dling of the observations, or of the way simulations are per- 
formed. In addition, a number of other observational evi- 
dence of MQ in nearby galaxies are now published, and we 
will discuss them in Section [6.31 



6.2 From a universal local star formation 

efficiency to global variations as a function of 
the host galaxy 

Assuming a universal conversion of dense gas into stars, our 
simulations reproduce variations of the global star forma- 
tion efficiency. Assuming that the small scale efficiency is 
universal seems reasonable based both on theoretical and 
observational grounds. A universal conversion of dense gas 
into stars is suggested by the observed linear correlation 
between SFR and mass of dense molecular gas (as mea- 
sured by the HCN luminosity) both in the Milky Way 
and other galaxies (|Gao SolomonI l2004l : IWu et al.|[2005bl : 
iLada et"aL 1201 j ). The global star formation efficiency then 
only depends on the fraction of molecular gas in a dense 
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phase, with for instance star burst ga laxies having higher 
HCN/CO ratios ('Gao Solo^^ l2004'V These observations 
can be interpreted by assuming a simple local law, consid- 
ering the dense clumps within giant molecular clouds as the 
standard units for star formation, with p roperties that do 
not change as a function of environment (IWu et alj l2005bl : 
iKrumholz ThomDSonll2007l : iKriimholz et al.ll2012l l. 

Instead, what changes from one galaxy to another is 
the amount of gas in these dense clumps, and their spa- 
tial distribution. For instance, in galaxy mergers, the g lobal 
efficiency of star fo rmation seems high er (IDaddi et alJ 2010: 
iGenzel et aDl2010l ). and simulations by Tevssie r et alJ (^20 IQ ) 
show that this can be explained by a universal local law but 
a gas density PDF skewed towards very high densities, with 
gas fragmented into very dense clumps. Var iations in the 
shape of the PDF have similarly been used bv lRenaud et al.l 
(|2012l ) to explain variations in global star formation effi- 
ciency from dwarf galaxies to high redshift mergers. 

Morphological quenching results from a similar effect: 
in ETGs the gravitational stability of gas discs is enhanced, 
and their density PDF is skewed towards lower densities. In 
some cases, they do not even reach the densities that would 
be typical of dense clumps (for instance gas that would be 
traced by HCN). In ETGs, the simulations thus predict a 
lower dense gas fraction than in spirals, which has tentatively 
been observed for a small sample in Paper XI. Given the 
shapes of the simulated PDFs, the HCN-to-Hi ratio should 
also be lower in ETGs than in spirals. 

If the picture of a universal local efficiency is true, then 
we also predict that ETGs should follow the same linear cor- 
relation between SFR and HCN luminosity as spiral galax- 
ies, which has been observed by iKrips et al.l (|2010h for 3 
SAURON ETGs but would require a larger sample to be defi- 
nitely confirmed. 



6.3 Other possible evidence of morphological 
quenching in the local Universe 

Numerous photometric studies, both at low and high red- 
shift, have shown that galaxies with lower levels of star for- 
mation tend to have a higher Sersic index, a higher stel- 
lar mass surface density (or 'inferred' velocity dispersion) 
and t o be more concentrated than star forming galaxies 
(e.g.. iKauffmann et a .1 l2003l: iFranx et al.1 [2008: Bel l 2008 ; 



IWilhams et al.1 l2010l : IWuvts et al.1 l201ll : iBeU et"aD [2012h . 
Combining the Fund amental Plane of ET Gs with stellar pop- 
ulation information, [Graves et al.l (|2009l l found that indeed 
the actual velocity dispersion is a better predictor a galaxy 
past star formation history than the galaxy stellar mass. The 
ATLAS project, thanks to accurate and unbiased masses 
from the dynamical modelling of the sample, confirmed that 
no other dynamical parameter pr edicts E TGs properties bet- 
ter than velocity dispersion (Cappellari 201 lb|, and Paper 
XX) . The same conclusion was reached by Wake et al.l 
using SDSS data. 

In Paper XX we used dynamical models to de-project 
the observed galaxy photometry and stellar kinematics. We 
demonstrated that velocity dispersion at fixed mass is a di- 
rect tracer of the bulge fraction in disc-like galaxies rather 
than the concentration of a spheroidal distribution, as fre- 
quently assumed. This suggests that the presence of a bulge 
is needed for quenching star formation, in agreement with 



the photometric study of iBell et al.1 (|2012l ). However, this 
does not mean that all these galaxies are undergoing mor- 
phological quenching; for instance the quenching could be 
due to the presence of an AGN (also correlated with bulge 
mass) . 

Only a study of the gas content of quenched galaxies can 
show if these galaxies are simply extremely gas poor, or in- 
stead contain gas but have a low SFR. Using the Arecibo 
radiotelescope, the atomic gas content of a large sample 
of low- redshift galaxies has recently be en measured by the 
GASS survey (|Schiminovich et al.ll2oTQh . Stacking of the HI 
spectr a for these ga l axies in bins of fixed stellar mass and 
color, iFabello e t al.l (}201l|) find that at fixed stellar mass 
and star formation rate or color, ETGs contain a lower 
mass of atomic gas than LTGs, a result which they claim to 
be in contradiction with morphological quench ing. However, 
the la ck of spatial resolution does not allow iFabello et al.l 
to determine the actual distribution of this atomic 
gas. Higher resolution observations show that, in ETGs, the 
HI is typically found in ve r y extended, low surface den sity 
discs (jMorganti et al.1 1 20061 : lOosterloo et alll2007l . lioTol . Pa- 
per XIII). In particular. Paper XIII shows that these galaxies 
lack the N(Hi)> 10^^ cm~^ atomic gas typical of the bright 
disc of spirals. Such low density gas would not be expected 
to form stars, independently of any morphological quenching 
effect. In addition, only measuring the atomic gas content of 
the galaxies is not enough to draw a conclusion on quench- 
ing, since Hi is generally ve ry poorly correlated with star 
formation (|Bigiel et al.ll2008l ). 

Molecular gas is much better correlated to SFR, as seen 
in the comparison of molecular gas mass and H/3 absorp- 
tion line strength for the ATLAS^^ sample (Paper IV). 
The COLD GASS survey (measuring CO mass for a sub- 
sample of the GASS galaxies) has reveale d strong trends 
between SFR. gas content and morphology ( Sainton ge et al.l 
l2011a| |b|,[20ll). They find that bulge-dominated galaxies fall 
into two categories. Some of them have very little cold gas 
and thus naturally do not form stars. The others contain gas 
(at surface densities similar to spiral galaxies) but tend to 
be inefficient at converting it into stars. More exactly, what 
(Saintonge et al. (2011b, 2012) find is a longer molecular gas 
depletion timescale (defined as Mhs/SFR) in galaxies with 
a high stellar mass surface density and a high concentration 
index, that is bulge- dominated galaxies. These galaxies are 
also clearly offset from the Kennicutt relation defined by spi- 
rals: the most bulge- dominated galaxies are all found below 
the mean relation. 



6.4 Implications for low and high redshift 

In agreement with others, our work shows that for a simi- 
lar gas content, ETGs form stars less efficiently than LTGs. 
This does not imply that red ETGs are gas-rich: in the Lo- 
cal Universe, only 22% of ETGs contain molecular gas, and 
when gas is pre sent the gas fract ion is often much lower 
than in spirals (I Young et al.|[201ll l. This is not in contra- 
diction with MQ, which is simply stating that when gas is 
present in ETGs it should be more stable against star for- 
mation. At high redshift [z > 1) galaxies are generally more 
gas-rich, l eaving more room for MQ to influence their SFR 
and color. ICeverino et "gH (|201 0') indeed found a case of MQ 
in one of their zoom cosmological simulations at z — 1.3, 
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Figure 8. Effect of a modified SFE on the g — r colour of sim- 
ulated galaxies. Following the method described in the text, we 
artificially reduce the SFE for a sample of simulated ETGs ex- 
tracted from a cosmological simulation that does not resolve disc 
stability and MQ. The black dashed line shows the initial distribu- 
tion of colours, with a lack of red galaxies. The red line shows the 
colour distribution when the SFE is reduced in ETGs as described 
in the text (the error bars show the range of values obtained for 
the different schemes we tested to reduce the SFE). Reducing the 
SFE by a factor of a few is enough to significantly redden the 
galaxies, and to create a red sequence. 



with a disc stabilized by the presence of a massive bulge. In 
their case, the bulge is formed by the coalescence of giant 
clumps during a phase of violent disc instability: this shows 
that MQ is an efficient mechanism whatever the origin of 
massive spheroids. 

However, this does not mean that all ETGs should 
be red and non star forming. Indeed, only relatively small 
amounts of gas can be kept stable. As already shown in M09 
and again in this paper, gas discs with too high surface den- 
sities can fragment, and then form stars with only a slightly 
lower efficiency than spiral galaxies. Another limitation is 
that for ETGs with a stable gas disc, interactions with small 
satellites could be enough to de-stabilize the disc and trigger 
some star formation. 

As a conclusion, MQ by itself cannot turn all ETGs 
red, but can definitely help in combination with other 
mechanisms. It also means that there might not be an 
absolute need t o prevent slow ha lo gas cooling through 
AG N feedback dCrotonejaD l2006l ) or gravitational heat- 
ing (| Johansson et al. 2009). These mechanisms may happen 
anyway, but in cases where they are inefficient and gas cools 
from the hot halo, this gas might simply join the disc and 
be stabilized, allowing the galaxy to remain red. 

To test the importance of MQ for the build-up of the 
red sequence, we use a hydrodynami c simulation of a cosmo- 
logical volume of (24/i"^ Mpc)^ (Gabor et al.ll201lh . This 
simulation captures galaxy accretion and outflows, but can- 
not resolve morphologies, gas stabilities, and SF efficien- 
cies. Producing red galaxies in that simulation required a 
"merger quenching" model in which an explicit superwind 
expels all the gas from the remnant immediately following 
a major merger. Because of gas re-accretion from the inter- 
galactic medium, the merger remnants still harboured sig- 
nificant levels of star formation and were not red enough. In 



add ition, no colour bim odality was found in the simulation 
(see iGabor et al.ll201lh . 

We aim to test if an explicit decrease of the SFE in 
these post-merger galaxies (i.e. imposing MQ in a simulation 
that is not able to resolve it) could produce more realistic 
galaxy colours and statistically change the properties of the 
sample. We select a small sample of 14 galaxies in the cos- 
mological simulation, targetting galaxies that have at z = 
a stellar mass between 10^° and 10^^ Mq and a relatively 
red colour (an intrinsicQ g — r greater than 0.45). We ex- 
tract their star formation histories and decrease their SFE 
whenever they have a bulge-dominated morphology. Unfor- 
tunately, the resolution of the simulation does not allow to 
follow morphological evolution with confidence, so instead 
we identify ETGs based on their g — r colour. We assume 
ETGs to have an intrinsic g — r greater than 0.5, which is 
the end of the blue cloud identified by Gabor et al. in that 
simulation, and roughly corresponds to galaxies with a spe- 
cific star formation rate (defined as the ratio between SFR 
and stellar mass) lower than 10~^^ yr~^. The galaxies bluer 
than 0.5 are unaffected by our SFE reduction scheme: they 
are assumed to be disc galaxies (some of them might be blue 
ETGs but not changing their SFE is then a conservative 
choice). The red galaxies are assumed to be ETGs, and have 
their SFE reduced by a factor going from 1 (if ^ — r = 0.5) to 
4 (if ^ — r > 0.65). 20% of the red galaxies are left unaffected 
to account for the fact that some of these galaxies could be 
discs. We also tested three other schemes for reducing the 
SFE in ETGs, but they all provide consistent results. 

As a function of time, as we reduce the SFE in some 
galaxies, we also modify their gas content self- consistently: 
a galaxy with a reduced SFE should consume less of its gas, 
thus providing additional fuel for future star formation. At a 
given time step, if we reduce the SFE, we then also increase 
the gas mass to keep the same total baryonic content as in 
the initial simulation. At the next time step, we use this 
new gas mass to compute a new SFR and g — r (the SFR is 
computed assuming the Kennicutt relation with the modified 
SFE). If the galaxy is still red, we continue to apply a SFE 
correction at the following time step, possibly by a different 
amount depending on the g — r colour. Note that we also 
track stellar mass loss: when the SFR is decreased, less gas 
is also re-injected into the ISM by stellar mass loss. 

We show in Figure [8] the resulting colour distributions 
for the original simulation and for the galaxies with a mod- 
ified SFE, using 7 snapshots between z — and z — 0.2 to 
increase the statistics (although all data points are then not 
totally independent). The error bars represent the range of 
values obtained with the four different schemes we tested. As 
expected, these schemes do not affect much the blue galaxies, 
but make the red ones significantly redder. In all cases, this 
creates a bimodal colour distribution. While this does not 
imply that MQ is the mechanism creating this bimodality 
in observed galaxies, and while our sample is small, it shows 
that a decrease by a factor of 2-5 of the SFE in ETGs is 
enough to have consequences for the reddening of galaxies. 



^ there is no dust correction 
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7 CONCLUSION 

Stability analyses predict that gas discs should be more sta- 
ble when they are embedded in an ETG instead of a spiral. 
This increased stability is due in part to the steeper potential 
well in ETGs, and mostly to the absence of stars contribut- 
ing to the disc self-gravity. We have previously shown in M09 
a cosmological zoom simulation in which a gas disc was nat- 
urally kept stable in an ETG, forming stars at a low rate, 
and allowing the galaxy to have red colors in spite of a large 
gas reservoir. 

In this paper, we further study the star formation effi- 
ciency in ETGs and spirals. Using both the simulation from 
M09 and a new sample of 5-pc resolution idealised simula- 
tions, we confirm that gas discs are significantly more stable 
in ETGs, and that their star formation efficiency is reduced 
by a factor of 2 - 5. From that small sample of simulations, it 
also seems that the slope of the Kennicutt relation is steeper 
for ETGs compared to spirals. This could be due to our star 
formation recipes, but is found both in the cosmological and 
idealized simulations, which are very different in terms of 
technique and resolution. Instead, it could also be a direct 
consequence of morphological quenching, with an even in- 
creased stability of gas disks at low densities, and at high 
densities a behaviour getting closer to what is found in spi- 
ral galaxies. A larger sample of simulations would be needed 
to understand that aspect. 

We also compare a sample of eigth local ETGs to twelve 
local spirals. For all these galaxies, we use a similar way 
of measuring CO and Hi surface densities, and SFR densi- 
ties from the non-stellar 8/im emission. We confirm that the 
ETGs are sighlty offset from the Kennicutt relation defined 
by the spirals, as also found for the COLD GASS sample by 
(jSaintonge et al.|[2012i ). We further validate our approach by 
performing a direct comparison between models and obser- 
vations. We run a simulation designed to mimic the stellar 
and gaseous properties of NGC524, a local lenticular galaxy, 
and find a gas disc structure and global star formation rate 
in good agreement with the observations. 

We thus conclude that given a similar average gas den- 
sity, the SFR is a few times lower in ETGs than in LTGs. 
This lower global efficiency arises in spite of a common uni- 
versal law for star formation on small scales. The cause is the 
different structure of the gas in both cases, with in ETGs a 
lack of gas at high densities. The difference between galaxy 
types is not only a difference of atomic-to-molecular gas con- 
version, but also (and most importantly) a difference within 
the molecular gas itself. Our simulations predict that in a 
quenched gas disc in an ETG the fraction of dense molecu- 
lar gas should be much lower than in spirals (with variations 
as a function of gas content and galaxy morphology). Such 
a lower fraction of dense molecular gas has tentatively been 
observed in Paper XL 
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